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Abstract —We consider the effects of decoding costs in energy 
harvesting communication systems. In our setting, receivers, in 
addition to transmitters, rely solely on energy harvested from 
nature, and need to spend some energy in order to decode 
their intended packets. We model the decoding energy as an 
increasing convex function of the rate of the incoming data. In this 
setting, in addition to the traditional energy causality constraints 
at the transmitters, we have the decoding causality constraints at 
the receivers, where energy spent by the receiver for decoding 
cannot exceed its harvested energy. We first consider the point- 
to-point single-nser problem where the goal is to maximize the 
total throughput by a given deadline subject to both energy 
and decoding causality constraints. We show that decoding costs 
at the receiver can be represented as generalized data arrivals 
at the transmitter, and thereby moving all system constraints 
to the transmitter side. Then, we consider several multi-user 
settings. We start with a two-hop network where the relay and the 
destination have decoding costs, and show that separable policies, 
where the transmitter’s throughput is maximized irrespective of 
the relay’s transmission energy profile, are optimal. Next, we 
consider the multiple access channel (MAC) and the broadcast 
channel (BC) where the transmitters and the receivers harvest 
energy from nature, and characterize the maxlmnm departure 
region. In all multi-user settings considered, we decompose our 
problems into inner and outer problems. We solve the inner 
problems by exploiting the structure of the particular model, 
and solve the outer problems by water-filling algorithms. 

Index Terms —Energy harvesting, throughpnt maximization, 
energy harvesting transmitters, energy harvesting receivers, de¬ 
coding costs, energy causality, decoding causality. 

I. Introduction 

Energy harvesting communications offer the promise of 
energy self-sufficient, energy self-sustaining operation for 
wireless networks with significantly prolonged lifetimes. En¬ 
ergy harvesting communications have been considered mostly 
for energy harvesting transmitters, e.g., m-m, with fewer 
works on energy harvesting receivers, e.g., lED-lMl- In this 
paper, we consider energy harvesting communications with 
both energy harvesting transmitters and receivers. 

The energy harvested at the transmitters is used for data 
transmission according to a rate-power relationship, which is 
concave, monotone increasing in powers. The energy harvested 
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at the receivers is used for decoding costs, which we assume 
to be convex, monotone increasing in the incoming rate 1^ . 
1321, lf35l - llMl . The transmission energy costs and receiver 
decoding costs could be comparable, especially in short- 
distance communications, where high rates can be achieved 
with relatively low powers, and the decoding power could be 
dominant; see 1^ and the references therein. 

We model the energy needed for decoding at the receivers 
via decoding causality constraints: the energy spent at the 
receiver for decoding cannot exceed the receiver’s harvested 
energy. We already have the energy causality constraints at the 
transmitter: the energy spent at the transmitter for transmitting 
data cannot exceed the transmitter’s harvested energy. There¬ 
fore, for a given transmitter-receiver pair, transmitter powers 
need now to adapt to both energy harvested at the transmitter 
and at the receiver; the transmitter must only use powers, and 
therefore rates, that can be handled/decoded by the receiver. 

The most closely related work to ours is ED, where the 
authors consider a general network with energy harvesting 
transmitters and receivers, and maximize a general utility 
function, subject to energy harvesting constraints at all ter¬ 
minals. Reference ETI carries the effects of decoding costs to 
the objective function. If the objective function is no longer 
concave after this operation, it uses time-sharing to concavify 
it, leading to a convex optimization problem, which it then 
solves by using a generalized water-filling algorithm. 

In this paper, we consider a similar problem with a specific 
utility function which is throughput, for specific network 
structures, with different decoding costs informed by network 
information theory. Eirst, we consider the single-user channel, 
and observe that the decoding costs at the receiver can be 
interpreted as a gate keeper at the front-end of the receiver that 
lets packets pass only if it has sufficient energy to decode. We 
show that we can carry this gate effect to the transmitter as a 
generalized data arrival constraint. Therefore, the setting with 
decoding costs at the receiver is equivalent to a setting with 
no decoding costs at the receiver, but with a (generalized) data 
arrival constraint at the transmitter ID. We also note that the 
energy harvesting component of the receiver can be separated 
as a virtual relay between the transmitter and the receiver; 
and again, the problem can be viewed as a setting with no 
decoding costs at the receiver but with a virtual relay with a 
(generalized) energy arrival constraint cni-iBi. 

We then consider several multi-user settings. We begin 
with a decode-and-forward two-hop network, where the relay 
and the receiver both have decoding costs. This gives rise 
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to decode-and-forward causality constraints at the relay in 
addition to decoding causality constraints at the receiver and 
energy causality constraints at the transmitter. We decompose 
the problem into inner and outer problems. In the inner 
problem, we fix the relay’s decoding power strategy, and 
show that separable policies are optimal ifTOl . Ifm . These 
are policies that maximize the throughput of the transmitter- 
relay link independent of maximizing the throughput of the 
relay-destination link. Thereby, we solve the inner problem 
as two single-user problems with decoding costs. In the outer 
problem, we find the best relay decoding strategy by a water¬ 
filling algorithm. 

Next, we consider a two-user multiple access channel 
(MAC) with energy harvesting transmitters and receiver, and 
maximize the departure region. We consider two different 
decoding schemes; simultaneous decoding, and successive 
cancellation decoding Each scheme has a different de¬ 
coding power consumption. For the simultaneous decoding 
scheme, we show that the boundary of the maximum de¬ 
parture region is achieved by solving a weighted sum rate 
maximization problem that can be decomposed into an inner 
and an outer problem. We solve the inner problem using the 
results of single-user fading problem 0. The outer problem is 
then solved using a water-filling algorithm. In the successive 
cancellation decoding scheme, our problem formulation is 
non-convex. We then use a successive convex approximation 
technique that converges to a local optimal solution iol, ED. 
The maximum departure region with successive cancellation 
decoding is larger than that with simultaneous decoding. 

Finally, we characterize the maximum departure region 
of a two-user degraded broadcast channel (BC) with energy 
harvesting transmitter and receivers. With the transmitter em¬ 
ploying superposition coding E2l . a corresponding decoding 
power consumption at the receivers is assumed. We again 
decompose the weighted sum rate maximization problem into 
an inner and outer problem. We show that the inner problem is 
equivalent to a classical single-user energy harvesting problem 
with a time-varying minimum power constraint, for which we 
present an algorithm. We solve the outer problem using a 
water-filling algorithm similar to the outer problems of the 
two-hop network and the MAC with simultaneous decoding. 

II. Single-User Channel 

As shown in Fig. [D we have a transmitter and a receiver, 
both relying on energy harvested from nature. The time is 
slotted, and at the beginning of time slot i £ {1,...,W}, 
energies arrive at a given node ready to be used in the same 
slot or saved in a battery to be used in future slots. Fet {Ei}fLi 
and {Ei}^i denote the energies harvested at each slot for 
the transmitter and the receiver, respectively, and let {pi}^i 
denote the transmitter’s powers. 

Without loss of generality, we assume that the time slot 
duration is normalized to one time unit. The physical layer is 
a Gaussian channel with zero-mean unit-variance noise. The 
objective is to maximize the total amount of data received 
and decoded by the receiver by the deadline N . Our setting 
is offline in the sense that all energy amounts are known prior 
to transmission. 
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Fig. 1. Single-user channel with an energy harvesting transmitter and an 
energy harvesting receiver. 


The receiver must be able to decode the kth packet by the 
end of the kth slot. A transmitter transmitting at power pi in 
the fth time slot will send at a rate g{pi) = ^ log 2 (1 +Pi), 
for which the receiver will spend (j){g{pi)) amount of power 
to decode, where f is generally an increasing convex function 
ED, ED, ES-ESl. In the sequel, we will also focus on 
the specific cases of linear and exponential functions, where 
(^(r) = ar + b, with a,b > 0, and ^{r) = c2‘^^ -b e, 
with c,d > 0 and c -|- e > 0. Continuing with a general 
convex increasing function we have the following decoding 
causality constraints for the receiver; 

k k 

'^Hg{Pi)) k = l,...,N (1) 

Therefore, the overall problem is formulated as; 

N 

max y2g{pi) 

p ^ 

1—1 

k k 

s.t. '^Pi<'^Ei, Vfc 

i—1 i=:l 

k k 

<Y. E„ Vfc (2) 

i—1 i—1 

where p denotes the vector of powers. Note that the problem 
above in general is not a convex optimization problem as ([D in 
general is a non-convex constraint since ^ is a convex function 
while 5 is a concave function E3l . Applying the change of 
variables g{pi) = and defining / = g~^ (note that f is a 
convex function), we have 

N 

max y Ti 

i=l 

k k 

s-t. '^f{n)<'^Ei, Vfc 

i—1 i—1 

k k 

Vfc (3) 

i—1 i—1 

which is now a convex optimization problem B3l . 

We note that the constraints in O- i-e-, < 

Sfci place upper bounds on the rates of the transmitter 
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by every slot k. This resembles the problem addressed in [T] 
with data packet arrivals during the communication session. 
In fact, when = r and Ei = bi, where bi is the amount 
of data arriving in slot i, these are exactly the data arrival 
constraints in m. A general convex (f) generalizes this data 
arrival constraint. We characterize the solution of Q in the 
following three lemmas and the theorem. The proofs rely on 
the convexity of / and (j) generalizing the proof ideas in ifTl . 

Lemma 1 {r*} is monotonically increasing. 




^n —1 


where 


In = arg mm <g 

Zti — I 




^n—1 






^ ^n—1 


with ig = 0, and n = 1,... ,JV. 


( 6 ) 


(7) 


Proof: Assume that there exists a time slot k such that > 
and consider a new policy obtained by replacing both 

rl and by fk = fk+i = ^ and observe that from 

the convexity of / and (j), we have 

f{h) + f{fk+i) < f{rl) + f{rl+i) ( 4 ) 

4){fk) + (l){fk+i) < + <^(rfc+i) (5) 

In addition, since both / and (j) are monotonically increasing, 
we have / {fk) < / (r^), and cj) {fk) < (p (r^). Therefore, the 
new policy is feasible, and can only save some energy either 
at the transmitter or at the receiver. This saved energy can be 
used to increase the rates in the upcoming time slots. Thus, 
the original policy cannot be optimal. ■ 

Lemma 2 In the optimal policy, whenever the rate changes 
in a time slot, at least one of the following events occur: 

1) the transmitter consumes all of its har\’ested energy in 

transmission, or 2) the receiver consumes all of its harvested 
energy in decoding, up to that time slot. 

Proof: Assume not, i.e., but both the transmitter 

and the receiver did not consume all their energies in the 
fcth time slot. Then, we can always increase and decrease 
without conflicting the energy causality or the decoding 
causality constraints. By the convexity of / and this 
modification would save some energy that can be used to 
increase the rates in the upcoming time slots. Therefore, the 
original policy cannot be optimal. ■ 

Lemma 3 In the optimal policy, by the end of the transmission 
period, at least one of the following events occur: 1) the 
transmitter’s total power consumption in transmission is equal 
to its total harvested energy, or 2) the receiver’s total power 
consumption in decoding is equal to its total harvested energy. 


Proof: Assume that both conditions are not met. Then, we can 
increase the rate in the last time slot until either the transmitter, 
or the receiver, consumes all of its energy. This is always 
feasible and strictly increases the rate. ■ 


Theorem 1 Let = 

the following 

r„ = min < g 


(j) A policy is optimal iff it satisfies 




^n—1 


Proof: First, we prove that the optimal policy satisfies (|6l) and 
©. We show this by contradiction. Let us assume that the 
optimal policy, that satisfies the necessary lemmas above, is 
not given by (|6l) and O and achieves a higher throughput. In 
particular, let us assume that it coincides with the policy given 
by (|6]l and (|2l) for all rates but has a different value 

for Let us denote the points of rate increase of this policy 
by {tfc}. Thus, there must exist a time index i' > in-i such 
that 


> min i g 




^ V - in-1 


^ ^n—1 


( 8 ) 


and let us consider two different cases. 

Assume that i' < in. If the transmitter’s energy is the 
bottleneck at i', then cannot be supported by the transmitter. 
On the other hand, if the receiver’s energy is the bottleneck 
at f, then cannot be supported by the receiver. Hence, ? „ 
is not feasible in both cases. Now, assume that i' > in. Then, 
there will exist a duration C + l,i'] where the rate has 
to decrease in order to satisfy feasibility. This violates the 
monotonicity property, and hence cannot be optimal. 

Second, let us show sufficiency. We show this again by 
contradiction. Let us assume that the policy that satisfies (|6l) 
and (I?]) is not optimal. In particular, let us assume that there 
exists another policy {r'} that coincides with it for all rates 
but has a different value for r„. Since this new policy 
should have higher throughput, we have r'n > rn- Now, assume 
i'n > in- Then, clearly is not feasible in the duration 
\in-i + l,*n]- On the other hand, if I'n < in, then by the 
monotonicity property, all upcoming rates {r'} for i > i'„ 
can only be larger than r^, which are all larger than r^. This 
makes the new policy infeasible by the end of slot in since 
consumes all feasible energy according to (|6]l and 0. Thus, 
the original policy is optimal. ■ 

Theorem [T| shows that decoding costs at the receiver are 
similar in effect to having a single-user channel with data 
arrivals during transmission and no decoding costs. This stems 
from the fact that the transmitter has to adapt its powers 
(and rates) in order to meet the decoding requirements at the 
receiver. Therefore, the receiver’s harvested energies and the 
function f control the amount of data the transmitter can send 
by any given point in time. 
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Fig. 2. Decoding costs viewed as a virtual relay. 


Alternatively, we can slightly change the single-user prob¬ 
lem Q by adding an extra variable fi as follows 

N 

max > fi 

r.r • ^ 

k k 

s.t. Vfc 

k k 

i—l i—1 

<ri, Vf (9) 

This gives the same solution as we will always have f* = r* 
satished for all i. Therefore, as shown in Fig. |2l we can view 
the single-user setting with an energy harvesting receiver, as 
a two-hop setting with a virtual relay between the transmitter 
and the receiver, with a non-energy harvesting receiver. To this 
end, we separate the decoding costs of the receiver, which are 
subject to energy harvesting constraints, as a relay which is 
subject to energy harvesting constraints in its transmissions, 
and consider the receiver as fully powered Uni-IBl. The 
receiver will only receive data if the relay has sufficient energy 
to forward it. In addition, this energy harvesting virtual relay 
has no data buffer, thus, its incoming data rate equals its 
outgoing data rate. The rate through this relay is controlled by 
Ei and </>. Thus, the decoding function cj) puts a generalized 
energy arrival effect to this virtual relay, in a similar way 
that it puts a generalized data arrival effect to the transmitter 
through Theorem [T] as shown in Fig. [T] 

It is worth mentioning that if we consider the special case 
where the receiver has no battery to store its energy, this will 
lead to the following decoding causality constraint 

^{g{p,))<E,, t = ( 10 ) 

which, in view of the generalized data arrival interpretation, 
can be modeled as a time-varying upper bound on the trans¬ 
mitter’s power in each slot 

pr = /(V'(^0) (11) 

where tpiEi) is the maximum transmission rate of a packet 
that Ei can handle at the decoder, and denotes its 

corresponding maximum transmit power. This problem has 
been considered in the general framework of m, and in Eol 


for the special case of a constant maximum power constraint. 
One solution for this problem is to apply a backward water- 
hlling algorithm that starts from the last slot backwards, where 
at each slot directional water-hlling 0 is applied only on 
slots whose maximum power constraint is not satished with 
equality. This might cause some wastage of water if the 
maximum power constraints are tighter than the transmitter’s 
energy causality constraints, which depends primarily on how 
the function (j) relates the transmitter’s and the receiver’s 
energies. 

III. Two-Hop Network 

We now consider a two-hop network consisting of a single 
source-destination pair communicating through a relay, as 
depicted in Fig. [3] The relay is full duplex, and it uses a 
decode-and-forward protocol. The relay has a data buffer to 
receive its incoming packets from the source. At the beginning 
of slot i, energies in the amounts of Ei, Ei, and Ei arrive at the 
source, relay, and destination, respectively. Unused energies 
can be saved in their respective batteries. 

Let ri and fi be the rates of the source and the relay, 
respectively, in slot i. Our goal is to maximize the total amount 
of data received and decoded at the destination by the deadline 
N. We impose decoding costs on both the relay and the 
destination. The problem is formulated as; 

N 

max A 

i=l 

k k 

S.t. (ri) <'^Ei, Mk 

i^l 

k k 

k k 

k k 

i=l i=l 

where the hrst constraint in (fT2l) is the source transmission 
energy causality constraint, the second one is the relay decode- 
and-forward causality constraint, the third one is the data 
causality constraint at the relay, and the last one is the 
destination decoding causality constraint. 
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Fig. 3. Two-hop energy harvesting system with both relay and destination decoding costs. 


We first note that that if the relay did not have a data 
buffer, the source and the relay rates will need to be equal, 
i.e., fi = Ti for all i. In this case, the problem reduces to be a 
problem only in terms of the source rates, and could be solved 
by straightforward generalization of the single-user result in 
Theorem [T] considering three constraints instead of two. In a 
sense, this would be equivalent to taking the effects of decode- 
and-forward causality at the relay and decoding causality at 
the receiver back to the source as two different generalized 
data arrival effects. This can be further extended to multi-hop 
networks with relays having no data buffers by taking their 
constraint effects all the way back to the source. 

In our setting, having a data buffer at the relay imposes non- 
obvious relationships among the source and the relay rates. To 
tackle this issue, we decompose the problem into inner and 
outer problems. In the inner problem, we solve for the source 
and relay rates after fixing a decoding power strategy for the 
relay node. By that we mean choosing the amounts of powers, 
the relay dedicates to decoding its incoming source 
packets. These amounts need to be feasible in the sense that 
< 'l2i=i Ei, V/c. This decomposes the decode-and- 
forward causality constraint into the following two constraints: 

k k k k 

E, - (5,, Vfc (13) 

i=l i=l i=l i=l 

In the next lemmas and theorem, we characterize the solution 
of the inner problem. The proofs of the lemmas are extensions 
of the ones presented in HD to the case of generalized data 
arrivals. 

Lemma 4 There exists an optimal increasing source rate 
policy for the inner problem. 

Proof: Assume that there exists a time slot k where rk > 
Tk+i- We have two cases to consider. First, assume fk > Pk+i- 
Let us define a new policy by replacing the fcth and fc -h 1st 
source and relay rates by r' = and f' = ^ 

respectively. By the convexity of / and and linearity of the 
data causality constraint, the new policy is feasible, and can 
only save some energy at the source or the relay. This energy 
can be used in later slots to achieve higher rates. 

Now, assume fk < ^k+i- We argue that the data arrival 
causality constraint is satished with strict inequality at time 


slot k. For if it were equality, we need to have fk > rk and 
fk+i < rk+i, which leads to rk < fk < fk+i < rk+i, an 
obvious contradiction. Now, we can find a small enough e > 0, 
such that defining a new policy by replacing the /cth and /c-|-lst 
source rates by — e and rk+i + e, respectively, we do not 
affect the relay rates. By the convexity of / and the new 
policy is feasible, and can only save some energy at the source. 
This energy can be used in later slots to send more data to the 
relay, and hence, possibly increasing the relay rates, and the 
end-to-end throughput. ■ 

Lemma 5 The optimal increasing source rate policy for the 
inner problem {r* } is given by the single-user problem solu¬ 
tion in (0 and m, where the transmitter’s and the receiver’s 
energies are given by {Ei} and {^i}, respectively. 

Proof: Let us denote the single-user solution by {rf}. Assume 
for contradiction that it is not optimal for the inner problem. In 
particular, let {r*} and {r'} be equal for i = 1,... ,k — l, and 
differ on the fcth slot. We again have two cases to consider. 
First, assume Ek > r'k. In this case, since by LemmalU {r *} is 
increasing, by similar arguments as in the proof of Theorem [D 
the policy {r* } will eventually not satisfy the source’s energy 
causality or the relay’s decoding causality constraints, at some 
time slot j > k. Hence, it cannot be optimal. 

Now, assume Ek < r}. We argue that this shrinks the 
feasible set of the relay’s rates. We show this by induction. 
By assumption of this case, it is true at time slot fc, that we 
have L* < Ei. Now, assume it is true that for some 

time slot j > k we have ^i, and consider 

the j -I- 1st time slot. If then we are back to 

the previous case where this cannot be feasible eventually. 
Therefore, the feasible set of the relay’s rates shrinks at time 
slot j -I- 1, and hence, shrinks all over k,... ,N. Thus, this 
case cannot be optimal either. ■ 

Lemma |5] states that the optimal source policy is sepa¬ 
rable Ha, HD in the sense that the source maximizes its 
throughput to the relay irrespective of how the relay spends 
its transmission energy. This stems from the fact that the 
relay has an infinite data buffer to store its incoming source 
packets. Therefore, once we fix a decoding power strategy at 
the relay, we get separability. The following theorem, which 
is an extended version of Theorem [D gives the optimal relay 
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rates for the inner problem. The proof is similar to that of 
Theorem [T] and is omitted for brevity. 


Theorem 2 Given the optimal source rates {r *}, the optimal 
relay rates for the inner problem is given by 


r„ = mm <g 




1 






i=l ^3 


'^n '^n—1 


'^n ^n—1 


(14) 

where in is the argmin of the expression in m as in (0-0, 
and io = 0. 


Denoting the solution of the inner problem by R{S), we 
now find the optimal relay decoding strategy {S *} by solving 
the following outer problem: 

max R{S) 

k k 

s.t. Vfc (15) 

i=l i=l 

We have the following lemma regarding the outer problem. 


Lemma 6 R{d) is a concave function. 


Proof: Consider two decoding power strategies ^ 2 , and 
let {ri,fi}, {r 2 ,f 2 } be their corresponding source and relay 
optimal inner problem rates, respectively. Let 6 g = 9Si + (1 — 
6 ) 62 , for some 0 < 0 < 1, and consider the rate policy defined 
by rg = Ori + (1 — 0 )r 2 , and rg = Ovi + (1 — 0 )r 2 , for the 
source, and the relay, respectively. By the convexity of / and 
(j), the policy {r^, f^} is feasible for the decoding strategy Sg. 
Therefore, we have 

N 

R{Sg) >J2r9^ = ORi^l) + (1 - 0 )R{S 2 ) (16) 

i=l 

proving the concavity of R{ 6 ). ■ 

Therefore, the outer problem is a convex optimization 
problem ll4^ . We propose a water-filling algorithm to solve 
the outer problem ini. We first note that R{S) does not 
possess any monotonicity properties in the feasible region. 
For instance, i?(E) = R{0) = 0, while R{S) is strictly 
positive for some 6 in between. Thus, at the optimal relay 
decoding power strategy, not all the relay’s decoding energy 
will be exhausted. To this end, we add an extra N + 1st slot 
where we can possibly discard some energy. We start by filling 
up each slot by its corresponding energy/water level and we 
leave the extra N + 1st slot initially empty. Meters are put 
in between bins to measure the amount of water passing. We 
let water flow to the right only if this increases the objective 
function. After each iteration, water can be called back if 
this increases the objective function. All the amount of water 
that is in the extra slot is eventually discarded, but may be 
called back also during the iterations. Since with each water 
flow the objective function monotonically increases, problem 



Fig. 4. Two-user MAC with energy harvesting transmitters and receiver. 


feasibility is maintained throughout the process, and due to 
the convexity of the problem, the algorithm converges to the 
optimal solution. 

IV. Multiple Access Channel 

We now consider a two-user Gaussian MAC as shown 
in Fig. m The two transmitters harvest energy in amounts 
and {E 2 i]f^i, respectively, and the receiver harvests 
energy in amounts receiver noise is with zero- 

mean and unit-variance. The capacity region for this channel 
is given by ll42l : 

ri < gipi) 

l 2 < gipT) 

ri + r2 < gipi + P2) (17) 

where pi and p 2 are the powers used by the first and the 
second transmitter, respectively. 

In addition to the usual energy harvesting causality con¬ 
straints on the transmitters 0, we impose a receiver decoding 
cost. We note that there can be different ways to impose this 
constraint depending on how the receiver employs the decod¬ 
ing procedure. In the next two sub-sections, we consider two 
kinds of decoding procedures, namely, simultaneous decoding, 
and successive decoding 1^ . Il42l . Changing the decoding 
model affects the optimal power allocation for both users so 
as to adapt to how the receiver spends its power. 

A. Simultaneous Decoding 

In this case, the two transmitters can only send at rates 
whose sum can be decoded at the receiver. A power control 
policy is feasible if the following are satisfied: 

k k 

'^Pii<'^Eii, Vfc 

k k 

'^P2i<'^E2i, Vfc 

i=l 2=1 
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k k 

'^4'{9{Pii+P2i)) <'^Ei, Vfc (18) 

i=l i=l 

From here on, we assume a specific structure for the decoding 
function (f) for mathematical tractability and ease of presenta¬ 
tion. In particular, we assume that it is exponential with param¬ 
eters c = 1, d = 2 and e = —1, i.e., (j){r) = g~^{r) = 2^’' — 1. 
Let Bj denote the total departed bits from the jth user by 
time slot TV. Our aim is to characterize the maximum departure 
region, V{N), which is the region of (i?i, B 2 ) the transmitters 
can depart by time slot TV, through a feasible policy. The 
following lemmas characterize this region 0. 

Lemma 7 The maximum departure region, 'D(N), is the 
union of all {Bi,B 2 ), over all feasible policies {pii,P 2 i}f=i, 
where for any fixed power policy, {Bi,B 2 ) satisfy 

N 

Bi < '^gipu) 

N 

B 2 < '^9{P2i) 

i=l 

N 

Bi + B 2 <'^g{pii+P2i) (19) 

i=l 

Lemma 8 I?(TV) is a convex region. 

Each point on the boundary of 'D{N), see Fig. |5] can be 
characterized by solving a weighted sum rate maximization 
problem subject to feasibility conditions (fTsT l. Let pi and p 2 
be the non-negative weights for the first and the second user 
rates, respectively. Assuming without loss of generality that 
Pi > p 2 , and defining p = we then need to solve the 

following optimization problem; 

N N 

X! 9{pu) + F X! sipii + P 2 i) 

2 = 1 2 = 1 


k k 

s.t. ^Pii<^Eii, yk 
2=1 2=1 

k k 

'^P 2 i<'^E 2 i, V/c 
2=1 2 = 1 

k k 

'^Pii+P 2 i <'^Ei, Vfc (20) 

2=1 2=1 

We note that the above problem resembles the one formulated 
in El for a diamond channel with energy cooperation. First, 
we state a necessary condition of optimality for the above 
problem. 

Lemma 9 In the optimal solution for by the end of the 
transmission period, at least one of the following occur: 1) 
both transmitters consume all of their harvested energies in 
transmission, or 2) the receiver consumes all of its harvested 
energy in decoding. 

Proof: Assume without loss of generality that transmitter 1 
does not consume all of its energies in transmission, and 
that the receiver also does not consume all of its energies in 
decoding. Then, we can always increase the value of piM until 
either transmitter 1 or the receiver consume their energies. This 
strictly increases the objective function. ■ 

We decompose the optimization problem (l20l i into two 
nested problems. First, we solve for p 2 in terms of pi, and 
then solve for pi. Let us define the following inner problem: 

N 

'^ 9 {pu+P 2 i) 

2=1 

k k 

'^P2z<'^Qz, Vfc (21) 
2=1 2=1 

where the modified energy levels Qi are defined as follows: 

Qi = Mi — Mi-i, 

Mi = min < ^ £' 2 ^ ,'^Ej -pij > , Mq = 0 (22) 

[i=i i=i J 

Then, we have the following lemma. 

Lemma 10 G(pi) is a decreasing concave function in pi. 

Proof: G is a decreasing function of pi since the feasible 
set shrinks with pi. To show concavity, let us choose two 
points p^^^ and Pi^\ and take their convex combination p® = 
0p^^^ -I- (1 — 0)Pi for some 0 < 0 < 1. Let p7^ and P 2 ^^ 
denote the solutions of the inner problem (ITTI) at p^^^ and Pi^\ 
respectively. Now, let p® = 6 *P 2 ^^ -I- (1 — 0)p^\ and observe 
that, from the linearity of the constraint set, p® is feasible with 
respect to pf. Therefore, we have 

G(p?) > '^g{Pu+P^i} 

2=1 


G(pi) = max 

P2 

S.t. 


max 

Pl:P2 











(25) 


N 

>^^9 {p[^i +P2i^) + (1 - &)9 [pu +P2i^) 

i^l 

= 6G (pi')) + (1 - e)G (pf)) (23) 

where the second inequality follows from the concavity of g. 

■ 

We observe that the inner problem (ISTT i is a single-user 
energy harvesting maximization problem with fading, whose 
solution is via directional water-filling of {Qi}fLi over the 
inverse of the fading levels {1 + Pii}fLi presented in [Jj. 
Next, we solve the outer problem given by: 

N 

max AiG(pi) 4-Vg(pH) 

Pi ^ 

1 = 1 

k k 

s.t. Vfc (24) 

i=l i=l 

where we define the water levels Ti = Li — Li_i, with Li = 
min ^*^1 E’j|, and Lq = 0. The minimum is 

added to ensure the feasibility of the inner problem. Note that, 
by Lemma [10] the outer problem is a convex optimization 
problem ||4^ . We first note that at the optimal policy, first 
user’s modified energies {Ti} need not be fully utilized by the 
end of transmission. This is because the objective function is 
not increasing in pi. To this end, we use the iterative water¬ 
filling algorithm for the outer problem proposed in Section HIJ 
to solve this outer problem. Since the problem is convex, 
iterations converge to the optimal solution. 

Note that the above formulation obtains the dotted points in 
the curved portion of the departure region in Fig. [5] Specific 
points in the departure region, e.g., points 1 and 3 in Fig. |3 
can be found by specific schemes ll45l . by solving the problem 
for the cases g,i = g ,2 and fiig ,2 = 0. 


B. Successive Cancelation Decoding 

We now let the receiver employ successive decoding, where 
it aims at decoding the corner points, and then uses time 
sharing if necessary to achieve the desired rate pair ll^ . Il42l . 
For instance, if the system is operating at its lower corner 
point, then the receiver first decodes the message of the second 
user, by treating the first user’s signal as noise, then decodes 
the message of the first user, after subtracting the second user’s 
signal from its received signal. For pi > ^ 2 , we are always 
at a lower corner point at every time slot, and therefore the 
weighted sum rate maximization problem can be formulated 
as: 

N ^ / 

max Ft y' 9{Pii) + ^2 V g ( , 
pi^P2 ^ ^ Vl-t-pi 

1=1 1=1 

k k 

s.t. '^Pu<'^Eii, Mk 

i=l i=l 

k k 

V/c 

i=l i=l 


+ Vfc 


where the last inequality comes from the fact that the receiver 
is decoding the second user’s message first by treating the 
first user’s signal as noise, and thereby spends (p (p (y^^)) 
amount of energy to decode this message, and then spends 
(p {g (pii)) amount of energy to decode the first user’s message 
after subtracting the second user’s signal. 

Observe that the last constraint, i.e., the decoding causality 
constraint, is non-convex. Therefore, one might need to invoke 
the time-sharing principle in order to fully characterize the 
boundary of the maximum departure region. In terms of the 
rates the problem can be written as: 


max 

ri,r2 


s.t. 


N 


N 


Pi ^ ru + M2 X! ^2* 

i=l i=l 

k k 

- Vfc 

i=l i=l 

k k 

2^’'“ (2^’'=^ -1) <'^E2i, Vfc 

i=l i=l 

k k 

-4 22’'=*-2 < Vfc 


i=l 


i=l 


(26) 


which is a non-convex problem due to the second user’s energy 
causality constraint. In fact, the above problem is a signomial 
program, a generalized form of a geometric program, where 
posynomials can have negative coefficients ll43l . Next, we use 
the idea of successive convex approximation ll40ll to provide 
an algorithm that converges to a local optimal solution. 

By applying the change of variables xji = 2^’'^* — 1, j = 
1 ,2, and some algebraic manipulations: 


N 


mm 

Xl,X2,tl,t2 


S.t. 


''2i 


i=l 

k k 

'^xii<'^Eii, Vfc 

i=l i=l 

k k 

'^{lEXii)X2i<^E2i, Vfc 

i=l i=l 

k k 

V/c 

i=l i=l 

tii ^ 1 “t“ Xii^ Vt 

t2z<i+X2i, Vt (27) 


Now, the problem looks very similar to a geometric program 
except for the last two sets of constraints. These constraints 
are written in the form of a monomial less than a posynomial, 
which will not allow us to write the problem in convex form 
by the usual geometric programming transformations m. 
We will follow an approach introduced in BTlI in order to 
iteratively approximate the posynomials on the right hand side 
by monomials, and thereby reaching a geometric program 
that can be efficiently solved ll43ll . Approximations should 
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be chosen carefully such that iterations converge to a local 
optimum solution of the original problem POl . Towards that, 
we use the arithmetic-geometric mean inequality to write: 

as, 


which holds for 0 < a < 1. In particular, equality holds at a 
point Xk > 0 if we choose a = Therefore, the monomial 

function u{x; ak) approximates the posynomial function 1 +x 
at X = Xk- Substituting this approximation, we obtain that at 
the k + lst iteration, we need to solve the following geometric 
program: 


min 

Xi,X2,ti,t2 


s.t. 


N 

i^l 

k k 

\/k 

i^l 

k k 

'^{1 + Xu)x 2 i <'^E 2 i, Vfc 

i=l 

k k 

'^Xu + X2i ^'^Ei, Vfc 

i=l i=l 



where 7 = 1,2, and x^^'^ is the solution of the 

fcth iteration. We pick an initial feasible point 
and run the iterations. The choice of the approximating 
monomial function u satisfies the conditions of convergence 
stated in iQi, and therefore, the iterative solution of problem 
(|29] | converges to a point (x^,X 2 ) that is local optimal for 
problem (l25l l. Finally, we get the original power allocations 
by substituting pl^ = x*,, and = (x*, -f 1) x^^. 


V. Broadcast Channel 

We now consider a two-user Gaussian BC with energy har¬ 
vesting transmitter and receivers as shown in Fig. |6] Energies 
arrive in amounts Ei, Eu, and E 2 i, at the transmitter, and the 
receivers 1 and 2, respectively. By superposition coding ll42l . 
the weaker user is required to decode its message while treat¬ 
ing the stronger user’s interference as noise. While the stronger 
user is required to decode both messages successively by first 
decoding the weaker user’s message, and then subtracting it 
to decode its own. The receiver noises have variances 1 and 
> 1 . 

Under a total transmit power P, the capacity region of the 
Gaussian BC is ll42l : 

?'1 < ^ log 2 (1 + aP), T 2 < i log 2 (^1 + (30) 

working on the boundary of the capacity region we have: 

P=[(t'^ - 1) 22’’2 -b =F (ri, ra) (31) 



Fig. 6. Two-user BC with energy harvesting transmitter and receivers. 


where F{r 1 ^x 2 ) is the minimum power needed by the trans¬ 
mitter to achieve rates ri and r 2 - Note that F is an increasing 
convex function of both rates. 

As in the MAC case, the goal here is to characterize the 
maximum departure region: 


max 

ri,r2 


s.t. 


N N 

fit ^ rii + Ai2 ^ r2i 
2=1 


k k 

'^F{ru,r2^) <'^Ei, 


2=1 2=1 

k k 


V/c 


+X2i) ^'^Eu, Vfc 
2=1 2=1 
k k 


^</>(r2*) Vfc 

2=1 2=1 


(32) 


where the first constraint in (l32l i is the source transmission 
energy causality constraint, and second and third constraints 
are the decoding causality constraints at the stronger and 
weaker receivers, respectively. Here also, we take the decoding 
cost function (j) to be (j>{r) = 2^’’ — 1. 

By virtue of superposition coding, we see that, in the 
optimization problem in (|3^ , the decoding causality constraint 
of the stronger user is a function of both rates intended for the 
two users, as it is required to decode both messages. While 
the decoding causality constraint for the weaker user is a 
function of its own rate only. By the convexity of F and </,, the 
maximum departure region is convex, and thus the weighted 
sum rate maximization in (l32l i is sufficient to characterize its 
boundary ||6l. In addition, the optimization problem in (l32l i is 
convex ll43l . 

We note that a related problem has been considered in ID, 
where the authors characterized transmission completion time 
minimization policies for a BC setting with data arrivals during 
transmission. There, the solution is found by sequentially solv¬ 
ing an equivalent energy consumption minimization problem 
until convergence. Their solution is primarily dependent on 
Newton’s method ll4Tl . Some structural insights are also pre- 
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sented about the optimal solution. In our setting, we consider 
the case with receiver side decoding costs, and generalize the 
data arrivals concept by considering the convex function (j). In 
addition, our formulation imposes further interactions between 
the strong and the weak user’s data, by allowing a constraint 
(strong user’s) that is put on the sum of both rates, instead of 
on individual rates. 


I'u ( = 0 , Vi 


22 ’-^^ - 1 - E 2 jj = 0 , Vi 
Vurii = 0, r]2ir2i = 0, Vi 


(38) 


We characterize the solution of the problem according to 
the relation between /ii and ^2 as follows. If fii > fi 2 , then 
due to the degradedness of the second user, it is optimal to put 
all power into the first user’s message. This way, the problem 
reduces to a single-user problem: 

N 

k k 

- 1 < Vfc (33) 

i=l i=l 

where the modified energy levels {Wi} are defined as follows: 

— Li Z/j — i, 

El, I , To = 0 (34) 


Lj = min 




i=i 


i=i 


max 

ri 


S.t. 


From here, we state the following lemmas 

Lemma 11 The sum rate is monotonically increas¬ 

ing. 

Proof: We prove this by contradiction. Assume that there 
exists some time slot k such that rik-\-r 2 k > ^i(fe-i-i)+^ 2 (fc+i)- 
From (|36] |, since the denominator cannot increase, the nu¬ 
merator has to decrease for the sum rate to decrease, i.e., 
Vik > VHk+i) > 0. From complementary slackness, we must 
have rik = 0. Therefore, in order for the sum rate to decrease 
we must have r 2 k > r 2 (fc+i), which in turn leads to p 2 fc = 0 . 

From (ITTI i. we know that for the weak user’s rate to 
decrease, the numerator has to decrease, i.e., we must have 
?? 2 (fc+i)-??i(fc+i) < mk-lik- Since r] 2 k = 0, this is equivalent 
to having ? 72 (fc+i) < pi(fc+i) — rjik- However, we know from 
above that rjik > ? 7 i(fc+i), i.e., ? 72 (fc+i) < 0 , an obvious 
contradiction by non-negativity of the Lagrange multipliers. 


On the other hand, if /ii < fi 2 , then we need to investigate 

the necessary KKT optimality conditions II43I . We write the Lemma 12 The weak users rate {r 2 j} is monotonically in- 
Lagrangian for the problem (l32l i as follows: creasing. 


N 


N 


-C = - /ri ^ rii - ^2 E ^2* 


2=1 


2=1 


N 


E^^ E(-^- 

i=i \j=i 


N 


2=1 \j = l 


N 


2=1 \i = l 


N 


N 


- E - E 


V 2 ir 2 i 


(35) 


Taking the derivative with respect to ru and r 2 i and equating 
to zero, we obtain: 


22 (rii-|-r 2 i) _ ~b T ?1 


^3 + 

M2 - Ml + 121 - Vli 


22 ^ 2 . ^ 


(36) 

(37) 


along with the complementary slackness conditions: 


A, I E “ 1) _ 0-2 - ] = 0, Vi 


Proof: We also prove this by contradiction. Assume that there 
exists some time slot k such that r 2 k > r 2 (fc+i). From 
(IJTI i. since the denominator cannot increase, the numerator 
has to decrease for the weak user’s rate to decrease, i.e., 
I? 2 (fc+ 1 ) - Mi(fe-I-I) < V 2 k - Vik- Let us consider two different 
cases. First, assume ijik > Mi(fe-i-i)- Therefore, we must 
have r] 2 k > ? 72 (fe+i) + (mu- - Vi{k+i)) > 0- and thus, 
by complementary slackness, r 2 k = 0, and hence, r 2 (k-ki) 
cannot be less since it cannot drop below zero. Now, assume 
Vik < Mi(fc-i-i)- In this case, by complementary slackness, 
ri(fc+i) = 0. By Lemma [m we have rik + r 2 k < r 2 (k+i), 
i.e., r 2 (fc+i) > r 2 k, which is a contradiction. ■ 

With the change of variables: pu = 2 ^Lii+»’ 2 i) _ 

P 2 i = — 1 , (l3^ becomes: 

N N 

max MiEff + (^ 2 -Mi)E5(M2*) 

Pt.P2 ^^ ^^ 

2=1 2=1 

k k 

s.t. ^(cr^ - l)p2i+Ptz < E^*’ 

2=1 2=1 

k k 

'^Pti<'^Eii, Vfc 
2=1 2=1 

k k 

E7'2*^E^2i, Vfc 
2=1 2=1 

Pti>P2i, Vi 


(39) 
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We now decompose the above problem into an inner and 
an outer problem and iterate between them until convergence. 
First, we fix the value of p 2 , and solve the following inner 
problem: 

N 

iJ(p 2 ) = max y^g{pti) 

pt ^ 

1—1 

k k 

s.t. '^pti<'^Vi, Vfe 

i=l 

Pti > P 2 z, Vi (40) 

where the modified energy levels are defined as follows 

^ = Bi — 

Bi = min \ - (tr^ - l)p 2 j > , Sq = 0 

J 

(41) 

We have the following lemma for this inner problem whose 
proof is similar to that of Lemma fTOl 


Lemma 13 iF(p 2 ) is a decreasing concave function in p 2 . 


We note that the P 2 vector serves as a minimum power 
constraint to the inner problem. Let us write the Lagrangian 
for the inner problem: 

N N / j j 

E = -Y^ g{pu) + X! S 

i—1 j — 1 \ 2=1 2 = 1 

N 

-^d‘i{Vti-P 2 i) (42) 

i=l 


Taking the derivative with respect to pu and equating to zero, 
we obtain: 


z2i=i^j Pi 


- 1 


(43) 


First, let us examine the necessary conditions for the opti¬ 
mal power to increase, i.e., pu < Pt{i+i)- This occurs iff 
Xi + /ii+i > Pi > 0. Thus, we must either have Ai > 0 
which means that, by the complementary slackness, we have 
to consume all the available energy by the end of the ith slot. 
Or, we have pi+i > 0 which means that Pt{i+i) = P 2 {i+i)- 
Next, let us examine the necessary conditions for the opti¬ 
mal power to decrease, i.e., pu > Pt(i+i)- This occurs iff 
Pi > Xi + Pi+i > 0 , and therefore, we must have pu = P 2 i- 
We note from Lemmas fTTI and fT 2 ] that both {p 2 i} and {pjj} 
are monotonically increasing. Therefore, we only focus on 
fixing an increasing feasible p 2 . This, when combined with 
the above conditions, leads to the following lemma. 


Lemma 14 For a fixed increasing p 2 , the optimal solution p( 
of the inner problem is also increasing. 


Proof: By the KKT conditions stated above, if we have 
Pti > Pt(i+i), then we must have pu = P 2 i- Thus, we will 


Algorithm 1 

1 : Initialize the status of each bin Si = Vi 

2 : Mark bins by their minimum power requirements {p 2 i}iLi 

3: Set k = N 
4: while /c > 1 do 

5: if Sk < P 2 k then 

6 : Pour water into the fcth bin from previous bins, in a 

backward manner, until equality holds 
7: else 

8 : Do directional water-filling over the current and 

upcoming bins {k, fc -|- 1,..., N} 

9: end if 

10 : Update the status of each bin 

11: fc ^ fc - 1 

12 : end while 


have Pt{i+i) < Pu = P 2 t < P 2 {i+i), i-e-, the minimum power 
constraint is not satisfied at the i + 1 st slot. ■ 

Therefore, choosing an increasing p 2 in the outer problem 
ensures that the inner problem’s solution p* is also increasing, 
and thereby, satisfies the conditions of Lemmas [TT] and [ 12 ] 
We solve the inner problem by Algorithm [T] The algorithm’s 
main idea is to equalize the powers as much as possible 
via directional water-filling 13 while satisfying the minimum 
power requirements. 

Observe that the algorithm gives a feasible power profile; 
it examines each slot, and does not move backwards unless 
the minimum power requirement is satisfied. If there is an 
excess energy above the minimum, say at slot fc, it performs 
directional water-filling which will occur if Sk > Sk+i (let 
us consider water-filling only over two bins for simplicity). 
Since the minimum power requirement vector p 2 is increasing, 
after equalizing the energies the updated status will satisfy 
Sk = Sk+i > P 2 {k+i) ^ P 2 k, i-e., the minimum power re¬ 
quirement is always satisfied if directional water-filling occurs. 
Also observe that the algorithm cannot give out a decreasing 
power profile since P 2 is increasing. 

According to the KKT conditions, the power increases from 
slot fe to slot fc-|-l only if ppk+i) = P 2 {k+i) or the total energy 
is consumed by slot fc. We see that the algorithm satisfies this 
condition. Power increases only if directional water-filling is 
not applied at slot fc, which means that either some of the 
water was poured forward in the previous iteration to satisfy 
Pt(k+i) = P 2 (k+i), or no water was poured which means that 
all energy is consumed by slot fc. 

A numerical example for a three-slot system is shown in 
Fig. |2| The minimum power requirements are shown by red 
dotted lines in each bin. According to the algorithm, we 
first initialize by pouring all the amounts of water in their 
corresponding bins. We begin by checking the last bin, and 
we see that it needs some extra water to satisfy its minimum 
power requirement. Thus, we pour water forward from the 
middle bin until the minimum power requirement of the last 
bin is satisfied with equality. This causes a deficiency in the 
middle bin, and therefore, we pour water forward from the 
first bin until the minimum power requirement of the middle 








12 


/ 


1: Initialization 2: Filling last bin first 

Fig. 7. Numerical example for the BC inner problem. 

bin is satisfied with equality. Since the problem is feasible, the 
amount of water remaining in the first bin should satisfy its 
minimum power requirement. In fact, in this example, there is 
an excess amount that is therefore used to equalize the water 
levels of the first two bins via directional water-filling. This 
ends the algorithm and gives the optimum power profile. 

We now find the optimum value of p 2 by solving the 
following outer problem: 

N 

^iH (P2) + ^g(P2*) 

k k 

Y.P2^<Y.K,, Vfc (44) 

i=l i=l 

where n = — ^ —, and the modified water levels iT, are given 

^ A‘2-/ri ‘ ® 

by: 

1 ; 

A, = min J ^ S 2 j , ^ ^ ^ j i , Ao = 0 (45) 

[j=l 3 = 1 3 = 1 J 

where the extra terms in the Ai expression are to ensure 
feasibility of the inner problem. By Lemma [13] the outer 
problem is a convex optimization problem HU. We solve it 
by an algorithm similar to that of the two-hop network outer 
problem, except that we only focus on choosing increasing 
power vectors p 2 in each iteration. By convexity of the 
problem, the iterations converge to the optimal solution. 

VI. Numerical Results 

In this section, we present numerical results for the con¬ 
sidered systems models. We focus on the specific case where 
g(x) = log(l + x), and </) = g~^. Starting with the single- 
user channel, we consider a five-slot system with energy 
amounts of E = [2,2,1,2.5,0.5] at the transmitter, and 
E = [1,1,0.5,2.5,3] at the receiver. The optimal rates 
in this case according to Theorem [1] are given by r* = 
[0.6061, 0.6061, 0.6061,1.2528,1.3863]. As we see, the rates 
are non-decreasing, which is consistent with Lemma [T] and 
they strictly increase only after consuming all the receiver’s 
energies in decoding by the end of the third slot, and again by 
the end of the fourth one, which is consistent with Lemma [2] 
In Fig. [8] we plot the maximum departure regions for a 
MAC with simultaneous decoding and successive cancellation 


3: Filling middle bin 4: Directional water-filling 

from first bin 


decoding. We consider a system of three time slots, during 
which the nodes harvest the energies: Ei = [0.5,1,2], 
E 2 = [1,2, 0.5], and E = [1.5, 2, 0.5]. We observe that the 
simultaneous decoding region lies strictly inside the successive 
decoding region. The latter, given by the geometric program¬ 
ming framework, is only a local optimal solution; one can 
therefore achieve even higher rates if a global optimal solution 
is attained. 

Finally, in Fig. [9] we provide some simulation results to 
illustrate the difference between the departure regions with 
and without decoding costs for a BC. We consider a system 
of three time slots, where the energy profile of the transmitter 
is given by E = [5,6, 7]. The maximum departure region with 
no decoding costs is shown in blue. We vary the energy profiles 
at the receivers to show the effect of the decoding costs on the 
maximum departure region. We start by setting Ei = [4, 5,6], 
and E 2 = [1,2,3], to get region A in red. Then we lower the 
values to Ei = [3,4, 5], and E 2 = [1,1.5, 2], to get region B 
in green. Finally, we lower the values again to Ei = [2, 3,4], 
and E 2 = [0.5,1,1.5], to get region C in brown. We note that 
as we lower the energy profiles at the receivers, the decoding 
causality constraints become more binding, and therefore, the 
region progressively shrinks. 

VII. Conclusions and Future Directions 

We considered decoding costs in energy harvesting commu¬ 
nication networks. In our settings, we assumed that receivers, 
in addition to transmitters, rely on energy harvested from 
nature. Receivers need to spend a decoding power that is 
a function of the incoming rate in order to receive their 
packets. This gave rise to the decoding causality constraints: 
receivers cannot spend energy in decoding prior to harvesting 
it. We first considered a single-user setting and maximized 
the throughput by a given deadline. Next, we considered two- 
hop networks and characterized the end-to-end throughput 
maximizing policies. Then, we considered two-user MAC and 
BC settings, with focus on exponential decoding functions, 
and characterized the maximum departure regions. In most of 
the models considered, we were able to move the receivers’ 
decoding costs effect back to the transmitters as generalized 
data arrivals', transmitters need to adapt their powers (and 
rates) not only to their own energies, but to their intended 
receivers’ energies as well. Such adaptation is governed by 
the characteristics of the decoding function. 


max 

P2 

s.t. 



































13 



Fig. 8. Departure regions of a MAC with simultaneous and successive 
cancellation decoding. 



Fig. 9. Depaiture regions of a BC with and without decoding costs. 

Throughout this paper, we only considered receiver de¬ 
coding costs in our models without considering transmit¬ 
ter processing costs. On the other hand, other works have 
considered the processing costs at the transmitter ll^ - ll2^ 
without considering decoding costs at the receiver. In their 
models, the transmitter spends a constant amount of power 
per unit time whenever it is communicating to account for 
circuitry processing; while in our model, the receiver spends 
a decoding power which is a function of the incoming data 
rate. As a future work, the two approaches can be combined 
to account for both the processing costs at the transmitter and 
the decoding costs at the receiver in a single setting. 
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